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Quantum entanglement shared between more than two particles is essential to foundational ques¬ 
tions in quantum mechanics, and upcor uing quantum information technologies. So far, up to 14 two- 
dimensional qubits have been entangledand an open question remains if one can also demonstrate 
entanglement of higher-dimensional discrete properties of more than two particles.^^ A promising 
route is the use of the photon orbital ang ular momentum (0AM), which enables implementa tion of 
novel quantum information protocolsand the study of fundamentally new quantum states.^^To 
date, only two of such multidimensional particles have been entangled^ albeit with ever increasing 
dimensionalityHere we use pulsed spontaneous parametric downconversion (SPDC)f^ to pro¬ 
duce photon quadruplets that are entangled in their QA M, or t ransverse-mode degrees of freedom; 
and witness genuine multipartite Dicke-type entanglement.^^^^ Apart from addressing foun dation al 
questions,^ this could find applications in quantum metrology, imaging, and secret sharing 


Twin photons that are created by SPDC are correlated 
in several degrees of freedom and exhibit quantum en¬ 
tanglement. Apart from the well-known polarization de¬ 
grees, the photons can also be correlated in their spa¬ 
tial degrees; this manifests itself in continuous wavevec- 
tor or (the Fourier-related) position entanglement.^^ We 
can also explore the spatial correlations using transverse 
paraxial optical modes, which can be chosen to form a 
discrete, orthogonal, and complete set, which is very use¬ 
ful for quantum information applications. In contrast 
to the 2 -dimensional polarization space, the transverse 
mode space is in principle infinite dimensional, only lim¬ 
ited by diffraction and the transverse size of optical com¬ 
ponents. Recently the use of computer-controlled holo¬ 
grams enabled steep increas e of th e single-particle Hilbert 
space up to 100 dimensionsAn experimentally use¬ 
ful choice of transverse modes are the Laguerre-Gauss 
(LG) modes where the azimuthal part factorizes and de¬ 
scribes phase vortices exp(i^ 0 ), where (j) is the azimuth 
and i = —OC...OC determines the twisting number of 
the wavefront; corresponding to an orbital angular mo¬ 
mentum of £h per photorP^ (in addition to the spin an¬ 
gular momentum). The LG and the related Hermite- 
Gauss modes have well-known propagation dynamics, 
thus they are suitable for long-distance distribution of 
high-dimensional entanglement. 

A photon pair produced by SPDG shows quantum en¬ 
tanglement in the transverse-mode and in particular the 
0AM degree of freedom,!^ this has already been used in 
a large number of quantum mechanical tests. The quan¬ 
tum correlations can be understood by considering the 
rotational symmetry of the SPDC process (on the wave¬ 
length scale, the crystal is rotationally symmetric); due 
to Noether’s theorem, this leads to conservation of 0AM. 
Therefore, we can write, in analogy to the well-known 
wave vector-space SPDC Hamiltonian, the 0AM SPDC 
Hamiltonian in case of a Gaussian pump beam as 
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H = -ihvh (1) 

£= — oo 

where aj is the creation operator for a photon with 0AM 
£ = —£ and k describes the strength of the nonlinear 
interaction. A single SPDC photon pair is produced by 
the lowest-order term of the series expansion of |T) = 
exp {—itIhH) | 0 ), leading to IT 2 ) = | 1 ^; 1 ^) which 

describes t wo p hotons that are perfectly anticorrelated in 
their OAIVP^ (l^/c) is a state with n photons in mode 
k). The single-pass amplitude gain j oc hzt depends lin¬ 
early on the pump beam intensity Ip. We explore here 
the next-order terms corresponding to the simultaneous 
production of 2 0AM photon pairs (for ^ 7 ^ 0 modes): 
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This state can be seen as a result of interference in a 
double-pair emission process, which for production of 
multiple photon pairs is much more stable than interfero¬ 
metric SPDCP^ experiments. Because the 4-photon term 
{Eq.§ depends quadratically on the pump beam inten¬ 
sity, we use a picosecond pulsed laser, further we use 1 nm 
wide band pass filters for the downconverted photons to 
limit spectral (and temporal) labelling.!^ Since we use 
type-I SPDC, all produced photons have the same polar¬ 
ization. We use nonpolarizing beamsplitters to separate 
them, and a combination of phase-only spatial light mod¬ 
ulation and projection onto the core of single mode fibers 
to perform projective measurements in transverse-mode 
space, see Fig. 

To study the SPDC produced state from Eq.[^ we record 
4-photon correlations for different detection modes. 
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Figure 1. Experimental implementation. A frequency- 
doubled (SFG) mode-locked picosecond laser is short-pass fil¬ 
tered (SPF) and focussed (LI) into the periodically poled 
potassium titanyl phosphate (PPKTP) crystal. The SPDC 
photons are spectrally filtered with a GaP plate and a band¬ 
pass filter (BPF), and distributed with beam splitters (BS) 
to the 4 equal detection units. The crystal facet is imaged 
with a telescope (L2 & L3) onto the spatial light modulators 
(SLM), which in turn is far-field imaged onto the core of the 
detection single mode fibers connected to single photon coun¬ 
ters. The pinhole (PH) selects the 1st diffraction order of the 
SLM holograms. We explore the 4-photon transverse-mode 
space by changing the holograms on the SLMs and recording 
4-fold coincidence events with a multi-channel time tagging 
computer card. 
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Figure 2. 4-photon OAM state production. Pump-power 
dependent single-photon count rate (a, b) and 4-fold coinci¬ 
dence rate (c, d). The single rates depend linearly on pump 
power, while the 4-fold coincidences rates show quadratic de¬ 
pendency. The comparison between zero time delay (a, c) and 
the case where photon 3 and 4 are delayed by 12.5 ns shows 
that most detected 4-photon detection events are indeed due 
to 4 photons that were produced within a single pump pulse. 



Figure 3. 4-photon OAM correlations. 4-photon joint de¬ 
tection probabilities in each of the 3 mutually unbiased bases: 
Hermite-Gauss {iL, H}, diagonal Hermite-Gauss {D, A}, and 
Laguerre Gauss {R, L}. Labels indicate the phase and am¬ 
plitude structure of the detected spatial modes. Gray bars: 
experimental data, error bars show statistical errors. Boxes: 
theory. Experimental integration time per data point was 24 
minutes. The probabilities are normalized to unity within 
each basis. 


pump powers, and relative detection times, shown in 
Fig.H We restrict ourselves to the fundamental Gaus¬ 
sian (G) mode and the 2-dimensional first-order mode 
space. The 3 mutually unbiased bases of the latter are 
{LG%,LGLi} = {R,L} in the Laguerre-Gauss ba¬ 
sis, {iLGo,i, = {H, V} in the Hermite-Gauss 

basis, and’ {HG^% HGi%} = {D, A} in the 45° ro¬ 
tated Hermite-Gauss basis, in analogy to the polariza¬ 
tion case.l^ We designed the experiment such that the 
single count rates depend only weakly on detection mode 
(Fig. I a, b). The 4-fold coincidence rate F for the case 
that all photons are produced by the same pump pulse 
(Fig. I c) is highest if all 4 photons are projected onto 
the fundamental Gaussian mode (|GGGG)), and depends 
quadratically on the pump power as expected. From the 
ratio between (d) and (c), whereas in (d) detectors (A, 
B) and (C, D) are set to detect photons produced in 
two different laser pulses, we can estimate the fraction 
of uncorrelated events to be 10% (F^^=^^-^^®/F^^^^). 
This is similar to the fraction of “forbidden” events (e.g., 
^HHHVHHHH ~ 0.1). Wc think the latter occur due 
to experimental imperfections; in contrast to polarization 
experiments we require here mode-matching between all 
4 detectors and the pump beam simultaneously. This 
argument is supported by the fact that both ratios are 
largely independent on pump power, suggesting that con¬ 
tributions from higher-order multiphoton states are low. 

These results support an intuitive explanation of the 
structure of Eq. keeping in mind that the photons are 
produced in pairs: The first term in Eq. [^that contains 
photons with different OAM \i\ occurs if the second pair 
is spontaneously emitted and uncorrelated to the first 
pair, while the second term corresponds to the case where 
the second pair is produced by a stimulated process giv¬ 
ing a perfect clone of the first pairThe different \i\ 
values in the first term in Eq. is then simply a conse¬ 
quence of forbidden perfect quantum cloning, both terms 
together demonstrate the possibility of optimal quantum 
cloning in stimulated SPDC.I^ 

Now we study the quantum correlations of the SPDC- 
produced 4-photon state, for which we focus on the 2- 
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dimensional first-order mode space with 0AM = ± 1 , 
and record 4-fold coincidences for all detector mode com¬ 
binations in each basis. Fig. [^compares experimental re¬ 
sults and theoretical prediction for the SPDC-produced 
state, which becomes in this case (Eq. |^: 


|^f)oc|2^;2,-) (3) 

with ^ = 1 . Note that Eq. is valid for any if we 
limit our detection to modes, a 2D Hilbert space per 
photon. Experimentally, the beamsplitters generate all 
possible permutations of photons (Eig. [^, so be- 

comes in the detector-basis 

\Df^) oc \e^E)^^M)+\ml)+\me)+\M^)+\Im) (4) 

which is the symmetric Dicke state of A" = 4 photons 
with A/2 = 2 excitations. This state is in particular 
interesting as it is robust to photon losses and has the 
largest distance from not genuine multipartite entangled 
states .1^ Entanglement in such states can easily be veri¬ 
fied by measuring the total collective spin along the x, y 
directionspl if = (Jj) + (J^) > 5 (for A = 4; 

Ji = \ where, e.g., = l(g)cr^(8)l(8)l), the 

state is non-separable, i.e., entangled. We obtain exper¬ 
imentally = 5.17 ± 0.09, thus verifying entangle¬ 

ment in the 4-photon 0AM state. 

/o') 

(W 4 ) can also be used to detect genuine multipartite 

entanglement if it violates + \/3 « 

5.23.1^ The proposed theoretical state in Eq.^ for which 
(IV 4 ) = 6 , is indeed genuine multipartite entangled, but 
our experimental result does not violate this bound. We 
argue that experimental imperfections are responsible: 
Apart from spectral-filtering issues, we have here the ex¬ 
treme requirement that all 4 detectors have to be mode 
matched simultaneously. In contrast to experiments on 
polarization entanglement, here, even small misalign¬ 
ment does not only reduce count rates but also alters 
the measurement projectors by inducing small rotations 
in the respective single-particle Hilbert space, and the 
4-fold mode-matching exponentially amplifies any mis¬ 
alignment. The experimentally obtained numerical val¬ 
ues should therefore be seen as a lower limit only. We 
can correct for this partially if we have access to the full 
density matrix, as we show now. 

We use tomographic quantum state reconstructiorP^ to 
obtain the most likely density matrix p describing the de¬ 
tected state. Experimental integration time is 120 s for 
each detector setting, and we additionally combine the 
data with the more accurate measurements of the basis- 
state correlations shown in Eig. In Eig. [^we compare 
the resulting experimental density matrix to the theoret¬ 
ical expectation (Eq. [^ . As a multipartite entanglement 
witness we use the one proposed in Ref. Ea which is also 



Figure 4. The reconstructed OAM quantum state. 

Modulus of the reconstructed experimental (a) and theoret¬ 
ical (b) density matrix (vertical axis scale: [0...0.3]). The 
discrepancy between experiment and theory is mainly caused 
by residual misalignment of the transverse-mode detectors. 


more resilient against noise compared to the collective- 
spin witness above. Violation of the inequality if [p] < 0 
indicates genuine multipartite entanglement of p. This 
witness is constructed such that it is optimal for the Dicke 
state with two excitations 1 ^ 4 ^^) and gives 1 in this case. 
For analysis of the experimental density matrix, we apply 
local single-qubit rotations and search for maximum vio¬ 
lation of the witness (see methods). For 3 independently 
measured data sets, we obtain [p] = {0.28, 0.32, 0.34}, 
thus giving strong indication of genuine 4-photon OAM 
entanglement. 

We have studied here the zero- and first-order modes 
with ^ = {0, ±1}. Exploration of the rich correlations 
in the full higher-dimensional multipartite state (Eq. 
will require higher pump beam intensities and therefore 
the use of nonlinear crystals with higher damage thresh¬ 
old, both within reach today. The 4-photon OAM en¬ 
tangled state that we have produced and characterized 
here might open new possibilities and protocols in multi¬ 
party quantum secret sharin^^^ with Dicke states in the 
sense that here, more information per photon quadruplet 
can be exchanged or the security increased by the high¬ 
dimensional nature of the OAM or transverse-mode de¬ 
grees of freedom.^ Eurther, the spatial correlations car¬ 
ried by our multi-photon states might enable new options 
in quantum metrology, microscopy, and imaging. 

Appendix 

Photon quadruplet creation in high-gain SPDC: 

We use a PPKTP crystal cut for degenerate collinear 
type-I SPDC. The photon pairs at 826 nm experience a 
different group index than the pump beam of 413 nm 
in PPKTP, therefore temporal labelling can occur. The 
group index difference is ng{2uj) — ng{uj) = 0.456 at 
300 K, corresponding to a group velocity dispersion of 
D = 1.5 ps/mm, and a group velocity walk-off length 
Lgy = 1.3 mm for At = 2 ps pulses. Thus, we choose 
a 1 mm long crystal. We separate the pump beam with 
an anti reflection coated GaP plate and select a narrow 
frequency window around the degenerate wavelength of 
826 nm with a 1 nm wide spectral filter. 
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The number of available entangled transverse modes is 
determined by the pump beam size in the crystal, at 
perfect phase matching. To achieve sufficient count rates, 
we focus the pump beam to 50 /im (while the detection 
mode waist in the crystal is 100 /im), which results in 
^ 9 transverse modes. 

If not noted otherwise, we use a pump power of 70 mW, 
which corresponds to a peak intensity at the focus of 
11.2 kW cm“^. This is close to the PPKTP damage 
threshold; above this, gray-tracking was observed. Ele¬ 
vated temperatures or other materials, such as period¬ 
ically poled Lithium Niobate, would enable the use of 
higher intensities. 

Entanglement witness optimization: 


n particles /^[p] optimized for a Dicke state with m 
excitations from Ref., 1 ^ only detect entanglement opti¬ 
mally in a particular basis, and optimization over basis 
rotations has to be done. Here, we optimize /^[p] by 
transformation of the experimentally determined (or the¬ 
oretical) density matrix p as follows: p ^ UpW^ where 
U = Ui 0 U 2 0 Us 0 U 4 with the generic unitary 


Ui = 


exp iai cos 7 ^ exp iPi sin 7 ^ 

-ex.pi{ai - 6i) sin-fi ex.pi{Pi - 6i)cosji 


for the real parameters Pi, Si, 7 ^}. These 16 param¬ 
eters are simultaneously optimized using robust uncon¬ 
strained numerical optimization routines. 


Local unitary transformations, complex basis rotations, 
are allowed to be applied to each part of a multipar¬ 
tite state without changing its entanglement properties. 
Most entanglement witnesses, including the witness for 
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